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We report a 75 As single crystal NMR investigation of LaFeAsO, the parent phase of a pnictide 
high T c superconductor. We demonstrate that spin dynamics develop a strong two-fold anisotropy 
within each orthorhombic domain below the tetragonal-orthorhombic structural phase transition at 
Tto ~ 156 K. This intermediate state with a dynamical breaking of the rotational symmetry freezes 
progressively into a spin density wave (SDW) below Tsdw ~ 142 K. Our findings are consistent 
with the presence of a spin nematic state below Tto with an incipient magnetic order. 

PACS numbers: 74.70.Xa, 76.60.-k 



The mechanism of high T c superconductivity in iron 
pnictides remains enigmatic [TJ. Earlier NMR measure- 
ments demonstrated that low frequency spin fluctuations 
associated with the spin density wave (SDW) instabil- 
ity grow toward T c near optimal doping [2HS], favoring 
the scenario of spin-fluctuation induced superconductiv- 
ity. There is, however, a complication; the slowing of 
the lattice vibrations accompanies that of spin fluctua- 
tions [6-8J. For example, neutron scattering measure- 
ments showed that LaFeAsO undergoes a tetragonal- 
orthorhombic structural phase transition at Tto ~ 156 
K, followed by a SDW ordering at T SDW ~ 142 K 
[6] [9j [10]. Moreover, the softening of the lattice begins 
at as high as ~ 200 K, and continues through Tto an d 
Tsdw down to Tuomains ~ 120 K [B], where the growth 
of the orthorhombic domains ends [LOl [TTJ . Theoretical 
analysis of LaFeAsO based on the frustrated J\- Ja model 
suggests that the Ising symmetry of Fe spins may be 
already broken below Tto without a three-dimensional 
magnetic long-range order 12J. Moreover, the inter- 
mediate temperature range between Tto an d Tsdw of 
LaFeAsO may be identified as a spin nematic state |13| . 

In such a nematic state, the spin correlations break 
the tetragonal symmetry, i.e. (S, • Si+ X ) = — (Si ■ Si+ y ) 
(we refer readers to Fig. 1 of a review article [T4J for a 
pictorial demonstration of the nematic state). More re- 
cent theoretical analysis based on an itinerant electron 
picture [Jjj] or an orbital fluctuation model [JB] also led 
to analogous conclusions. The prospect of observing such 
a magnetic analogue of a liquid crystal below Tto with 
an incipient ("fluctuating") magnetic order, and its po- 
tential link with the mechanism of high T c superconduc- 
tivity, has stimulated strong interest among researchers. 
The past experimental efforts searching for the signature 
of nematicity were focused primarily on the BaFe2As2 se- 
ries (e.g. [17H20] ). However, the proximity between the 
structural and SDW transitions, and/or the twinning of 
orthorhombic domains hampered these efforts. 



In this Letter, we report a microscopic 75 As NMR in- 
vestigation of LaFeAsO for a single crystal |UJ, and com- 
pare our results with neutron scattering [10] and mag- 
netic susceptibility x measured for the same piece of ~ 20 
mg crystal. The usage of a single crystal enabled us to 
resolve complicated changes of NMR lineshapes across 
Tto and Tsdw for the first time, and find the signature 
of the spontaneous breaking of the rotational symme- 
try. We will demonstrate that low frequency spin dy- 
namics indeed exhibit a strong anisotropy within each 
orthorhombic domain with a two-fold symmetry below 
Tto ~ 156 K. Moreover, the anisotropic spin state freezes 
progressively into a static SDW from Tsdw ~ 142 K to 
^Domains ~ 120 K, and the SDW ordered and paramag- 
netic domains coexist in a broad range of temperature. 
Our findings uncover the presence of an unconventional 
intermediate spin state below Tto with the signatures of 
spin nematicity. 

In Fig. 1, we summarize representative field-swept 75 As 
NMR lineshapes observed at Jnmr — 58.159 MHz (nu- 
clear spin / = 3/2). When we apply an external mag- 
netic field B ext along the c-axis, a sharp paramagnetic 
(PM) central peak appears at B^f er ~ 7.96 T for the 
I z = +1/2 to —1/2 transition, as shown in Fig. 1(b); the 
resonant condition is Jnmr = (1+ 75 K)-f n B^ ter , where 
the nuclear gyromagnetic ratio 7„/27r = 7.2919 MHz/T. 
75 K ~ 0.002 (~ 0.2 %) is the Knight shift, which mea- 
sures the product between the local spin susceptibility 
and the hyperfine coupling constant. We note that the 
full-width at half-maximum of the PM central peak is 
as sharp as 6.5 kHz (8.9 Oe) at 290 K in the absence of 
orthorhombic distortion. The narrow linewidth is com- 
parable to that of undoped BaFe2As2 [351 US], an d attests 
to the high homogeneity of our crystal. 

In Fig. 1(a), we observe an I z = ±3/2 to ±1/2 param- 
agnetic satellite peak near Bsateihte _ g j ^ (connected 
by a dotted line). We also found the strongly tempera- 
ture dependent antiferromagnetic (AF) central peak aris- 
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FIG. 1: (Color online) Representative 75 As NMR lineshapes 
observed at 58.159 MHz. (a) The paramagnetic (PM) satellite 
peak (o), and the antiferromagnetic (AF) central peak marked 
with a downward arrow (•) in B ext \\ c. For clarity, the verti- 
cal axis is shifted for different temperatures. The dotted line 
indicates the shift of the PM satellite peak with temperature. 
The horizontal axis at the top of (a) measures the split (-B£/) 
from the PM central peak at 7.96 T in (b). (b) The PM cen- 
tral peak with B ext || c (upper right) and B ext \\ ab (lower 
left). For clarity, the origin of the vertical axis is shifted for 
B ext || c. The formation of twinned orthorhombic domains 
along the two orthogonal axes, as schematically shown in Fig. 
2(a), results in the splitting of the B ext \\ ab lineshape below 
T TO ~ 156 K. 



ing from the statically SDW ordered domains (marked by 
downward arrows) , but only below the onset temperature 
of Tggffi » 135 K (< T SDW ). The AF central peak is 
shifted from the PM central peak at Bf x f er ~ 7.96 T 
in Fig. 1(b) by a static hyperfine field along the c-axis, 
{B%f). (B^f) originates from the ordered Fe magnetic 
moments M Fe , and {B c hf ) oc M Fe [22]. In Fig. 2(b), we 
deduce the temperature dependence of (B^f) from Fig. 
1(a), and compare the results with the neutron scattering 
data of M Fe [10]. 

The splitting between the central and satellite peaks, 
Vq = in{Bf x f er - B s e f t elllte ) » 9.2 MHz, measures the 
nuclear quadrupole interaction with the electric field gra- 
dient (EFG). We summarize the temperature dependence 
of Vq in Fig. 2(c). The EFG is the second derivative of 
the Coulomb potential arising from electrons and ions 
near the observed 75 As sites. Below Tto ~ 156 K, v"q 
exhibits a sharp downturn, because the EFG is sensitive 
to the local structural environment. Earlier diffraction 
measurements showed that the growth of orthorhombic 
domains finally comes to an end at Tuomains w 120 K 
[TTJ1 fTT] ; vq also levels off below Trj omains ■ 

In Fig. 1(b), we also show representative NMR line- 
shapes with B ext applied along the tetragonal a-axis 
within the FeAs planes. Above Tto, an tetragonal do- 
mains are equivalent, hence we observe a single central 



FIG. 2: (Color online) (a) Schematics of the orthorhombic 
domain with B ex t \\ b (upper left) and B ex t || a (upper right) 
within twinned ab-planes. Bottom: rotation by azimuthal 
angle 0. (b) (•): (B c hf ) deduced from Fig. 1(a). (A): the 
sub-lattice magnetization Mp e deduced as the square-root of 
the magnetic Bragg scattering intensity (after ref. 10 ), nor- 
malized at 4.2 K. (c) Vq measured for the paramagnetic (•) 
and antiferromagnetic (■) peaks, (d) The magnetic suscep- 
tibility Xat (•) and Xc (A) measured in 1 T with dc-SQUID. 
Also shown is the second derivative of \ a b (□)• Notice the 
presence of two kinks in \ab at ~ 156 K and ~ 133 K, as 
evidenced by the negative maxima of d 2 Xab/dT 2 . 



peak. Below Tto, the elongated a-axis of each domain 
points either along the direction of B ext , or orthogonal to 
B extl due to twinning of the orthorhombic domains, as 
sketched in Fig. 2(a). Moreover, orthorhombic distortion 
breaks the axial symmetry of the EFG. The difference in 
the second order effects of the nuclear quadrupole inter- 
action thus results in splitting of the central peak; the 
peak positions depend on the direction of the elongated 
a-axis relative to B ext . We confirmed that the NMR 
line splitting in Fig. 1(b) is caused entirely by the dif- 
ference in the second order nuclear quadrupole effects, 
which is inversely proportional to B ext [24J. That is, 
the NMR Knight shift is still axially symmetric, and 
shows very little temperature dependence below Tto, 
75 K a = 75 K b ~ 0.22 ± 0.03 %. In view of the fact 
that expansion of the lattice from Tto to 290 K results 
in a larger value of the quadrupole frequency as shown in 
Fig. 2(c), we tentatively assign the peak near B ext ~ 7.91 
T to the orthorhombic domains with B ext | a; the other 
peak near B ext ~ 7.92 T arises from the orthorhombic 
domains with B ext \\ b instead. From the splitting of an 
I z = ±3/2 to ±1/2 satellite peak, we estimate Vq ~ 5.2 
MHz and v b Q ~ 3.9 MHz at 145 K. We also confirmed 
that these double peaks collapse into one when we ap- 
ply B ext along the [110] direction within the ab-plane, 
because B ext points along the diagonal direction for all 
orthorhombic domains in such a geometry. 

For bulk averaged measurement techniques such as re- 
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FIG. 3: (Color online) (a) Temperature dependence of 1/T±T 
measured for the PM central peak with B ext along the [110] 
(•), and the c-axis (A) direction. The solid curves are guides 
to the eye only, (b) \jT\T as a function of the field orien- 
tation ij> as defined in Fig. 2(a), within each orthorhombic 
domain. The results are symmetrized for <f> < due to crystal 
symmetry. Below Tto ~ 156 K, 1/TiT exhibits increasingly 
strong dependence on the in-plane orientation of B ext , which 
can be modeled by the function 1/TiT = A + Bcos 2 <f> (solid 
curves). The accidental superposition between different NMR 
peaks makes accurate measurements of 1/TiT unfeasible be- 
low 135-140 K, depending on the direction of B ext - 

sistivity, the twinning of orthorhombic domains would 
result in experimental data averaged over two orthogo- 
nal directions within the ab-plane, unless one applies a 
uniaxial stress [T7H2U]. It is not straightforward, then, 
to probe the spontaneous breaking of the rotational sym- 
metry. In contrast, NMR is a local probe. Since we 
have succeeded in resolving the central peaks in Fig. 1(b) 
for different orientations, we can investigate the in-plane 
anisotropy of spins within each orthorhombic domain. 

In Fig. 3(a), we summarize the temperature depen- 
dence of the nuclear spin-lattice relaxation rate 1/Ti 
divided by T, l/T x T, for the PM central peaks [25] . 
Our preliminary \jT\T results for the SDW ordered 
state are very similar to the case of BaFe2As2 [22], 
and beyond the scope of the present work. \jT\T 
measures low-frequency spin fluctuations: (1/T 1 T) Q cx 
£ q |F(q)| 2 x"(q, Inmr)/ Inmr., where a specifies the di- 
rection of B ext , |f(q)| 2 and x" are the hyperfine form fac- 
tor [26] and the imaginary part of the dynamical electron 
spin susceptibility, respectively, and the q summation is 
taken within the Brillouin zone. \fT\T shows a mild in- 
crease with decreasing temperature down to Tto due to 
the slow growth of a short range SDW order. Once we en- 
ter the orthorhombic phase below T^Oi 1/T\T measured 
with B ext || [HO], a-, and b-axis begins to show rapid 
growth. This implies that the orthorhombic distortion 
enhances low frequency components of antifcrromagnetic 
spin fluctuations, and that the dynamic SDW is rapidly 
slowing down [31 [231 121] • This conclusion is consistent 



with the downturn of x below Tto [S], also observed for 
our crystal as shown in Fig. 2(d). At first glance, the 
enhancement of \jT{T below Tto is much weaker with 
B ext || c. This is simply because the transferred hyper- 
fine fields at 75 As sites from their four nearest-neighbor 
Fe sites are geometrically cancelled out within the ab- 
plane in this configuration, i.e. |F(Q)| 2 = for the SDW 
ordering wave vectors Q, and the contributions of AF 
spin fluctuations to (l/T{T) c are "filtered out" [25]. The 
growing anisotropy of \/T\T between the ab- and c-axis 
orientations observed below Tto therefore has little to do 
with that of the critical dynamics of Fe spins near Tsdw- 
Next, let's turn attention to the angular dependence of 
\jT\T within the ab-plane, which has been proposed as 
a novel probe of spin nematicity ]26] . As summarized in 
Fig. 3(b), we don't observe any ^-dependence of \jT\T 
above Tto- Once we enter the orthorhombic phase below 
Tto i 1/TiT begins to develop a strong anisotropy within 
each orthorhombic domain. The anisotropy reaches as 
much as a factor of ~ 2 by < 140 K. In view of the 
very small difference between the lattice constants of the 
a- and b-axis (~ 0.5 %) [91 [U3 IS], our finding is quite 
unexpected for paramagnetic spin fluctuations. 

As explained above, the Knight shift remains axially 
symmetric within experimental uncertainties below Tto ■ 
It is therefore unlikely that the uniform spin suscepti- 
bility or the hyperfine form factor |f(q)| 2 develops a 
sizable anisotropy within the ab-plane below Tto- We 
therefore conclude that low frequency Fe spin dynam- 
ics as reflected in x" (cL Jnmr) locally develop a strong 
rotational anisotropy by a factor of ~ 2 within each 
orthorhombic domain below Tto, without exhibiting a 
three-dimensional magnetic order. We note that the in- 
tensities of the a- and b-peaks in Fig. 1(b) are compa- 
rable, hence FeAs planes are randomly twinned. When 
averaged over the entire single crystal, Fe spin dynamics 
would appear almost isotropic within the ab-plane. 

The anomalous behavior of Fe spins below Tto is not 
limited to the in-plane anisotropy of their dynamics. In 
Fig. 4, we summarize the temperature dependence of the 
integrated intensities of the NMR signal. An unusual 
aspect of the signal intensity is that the paramagnetic 
NMR peaks don't disappear suddenly at Tsdw ~ 142 
K, and linger well below Tsdw down to T Domains ~ 120 
K. On the other hand, AF NMR signals from statically 
SDW ordered domains emerge progressively, only below 
t sdw ~ 135 K (< t sdw)- Notice that the PM central 
peak is still clearly observable below Tsdw in Fig- 1(b). 
Our finding that the PM signal intensity is as large as 
- 60 % at Tgp'l^ implies that ~ 60 % of the sample 
volume remains paramagnetic at 135 K. That is, AF and 
PM domains coexist even below Tg-pffi, although neu- 
tron scattering begins to detect magnetic Bragg peaks 
below Tsdw [10] ■ In addition, 1/TiT measured for the 
residual PM peak does not blow up at Tsdw, an d con- 
tinues to increase. If the SDW ordering in our LaFeAsO 
single crystal was a typical second order phase transi- 
tion, \jT\T would diverge at Tsdw due to the critical 
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FIG. 4: (Color online) The temperature dependence of the 
PM and static AF volume fractions at the NMR time scale, 
measured as the integrated intensity of the PM and AF central 
peaks (normalized by the Boltzmann factor). The solid curves 
are guides to the eye. Notice that the PM phase coexists 
with the AF phase even below Tsdw down to Toomains- We 
measured the spin echo with a delay time 2r = 40 fis; the 
transverse relaxation time T2 for the observable NMR signals 
is always long (~msec), and does not appreciably affect the 
intensity. Notice that lowering B ex t from ~ 7.96 T to ~ 3.85 
T does not alter the qualitative features. 



slowing down of spin fluctuations, followed by a sudden 
disappearance of PM NMR signals below Tsdw- 

The reason NMR and neutron scattering detect differ- 
ent onset temperatures of the SDW is that each experi- 
mental probe has a different characteristic measurement 
time scale [251 US] ■ Elastic neutron scattering measure- 
ments would consider the SDW 'static' when fluctuations 
slow down to below the instrument resolution of ~ lmeV. 
This means that neutron scattering can take an instanta- 
neous picture of Fe spins with a 'shutter speed' of ~ 10~ n 
sec even if they are still slowly fluctuating. In contrast, 
AF NMR signals in Fig. 1(a) become observable only 
when Fe spins in the SDW become static, to the ex- 
tent that (Bfrf) is time-independent over the duration 
of our spin echo measurements, 40 /isec. If Fe moments 
are fluctuating faster than 40 /xsec in some segments of 
FeAs planes, NMR would see them as motionally aver- 
aged out (i.e. paramagnetic). In other words, the dif- 
ferent onset temperatures of the SDW ordering between 
neutron and NMR data indicate that the fluctuations of 
the SDW continue to slow down from Tsdw ~ 142 K to 
t sdw ~ 135 K - The coexistence of PM and AF domains 
below 142 K implies that the fluctuation time scales have 
a broad distribution throughout the FeAs planes. We 
also found that Xab exhibits an additional kink at ~ 133 
K, as shown in Fig. 2(d). Since SQUID measures the 
time-independent response of spins, the acceleration in 
the suppression of Xab below Tgp ff is consistent with 
NMR. 

It is interesting to realize that these unusual behav- 



iors of LaFeAsO share similarities with the spin stripes 
in high T c cuprates [30] , where the concept of namaticity 
was originally proposed for unconventional superconduc- 
tors |31j . In the striped cuprates, the spin stripes pro- 
gressively slow down below a charge ordering at ~ 70 
K [30] (instead of Tto)\ elastic neutron scattering, ^SR, 
and NMR detect the emergence of a static SDW at their 
respective measurement time scale below ~ 50 K |30) . 
~ 30 K [32, and ~ 1.6 K [2E], respectively. Moreover, 
PM NMR signals linger well below - 50 K [28] . 

The temperature dependence of the AF NMR signal 
intensity in Fig. 4 indicates that the volume fraction of 
the static SDW at the time scale of NMR gradually in- 
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K. Once the growth of the orthorhombic domains ends 
at TDomains [El HI], the intensity saturates at ~ 60% 
of the paramagnetic intensity above Tsdw- The missing 
signal intensity suggests that (B^t) is still modulating 
and/or T2 is too fast for a spin echo NMR signal to form 
in ~ 40 % of the sample volume. Our finding is consistent 
with the earlier /iSR measurements for a powder sample 
of LaFeAsO; the muon precession with a well-defined fre- 
quency of ~ 23 MHz takes place only in ~ 70 % of the 
sample volume, and ~ 30 % of the [iSR signal is strongly 
damped [33]. It remains to be seen whether these ul- 
tra slow dynamics of the SDW in 30-40 % of the sample 
volume are caused by the motion of the anti-phase do- 
main boundaries [5^ and/or the finite size effects of the 
orthorhombic domains. 

To summarize, we demonstrated that Fe spin fluctua- 
tions in a LaFeAsO single crystal begin to slow down be- 
low Tto ~ 156 K, accompanied by the local breakdown 
of the rotational symmetry of spin fluctuations in each 
of the randomly twinned orthorhombic domains. The 
way that the static SDW develops is also unconventional. 
The paramagnetic and SDW ordered domains coexist in 
a wide range of temperature below Tsdw ~ 142 K due to 
a distribution in the fluctuation time scales of the SDW. 
A large volume of FeAs planes sees freezing of the static 
SDW only below TDomains ~ 120 K. Our findings point 
towards the presence of a novel spin state between T T q 
and TDomains with a dynamically broken rotational sym- 
metry and an incipient magnetic order, i.e. the elusive 
spin nematic sate. 
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